Time delay matrix at the spectrum edge and the minimal chaotic cavities 
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Using the concept of minimal chaotic cavities, we give the distribution of the proper delay times 
of Q = —ikS" 1 §J§ at the spectrum edge with a scattering matrix S belonging to circular ensembles 
CE. The three classes of symmetry (/? = 1,2 and 4) will be analyzed to show how it differs from the 
distribution obtained in the bulk of the spectrum. In this new class of universality at the spectrum 
edge, more attention will be given to the Wigner's time t w = Tr(Q) and its distribution will be 
given analytically in the case of 2 modes scattering. The results will be presented exactly at all 
the Fermi energies without any approximation. All this will be tested numerically with an excellent 
precision. 



PACS numbers: 72.10.-d 05.45.-a, 73.23.-b, 



Since it was introduced by Eiscnbud and Wignerfl], the 
time delay in scattering problems was extensively studied 
and used in describing quantum transport. This quan- 
tity was generalized by Smith to the notion of time delay 
matrix Q = —ihS^^ in multichannel scattering[2]. The 
eigenvalues of the matrix Q, known as the proper delay 
times play an important role in considering interaction 
in low dimensional systems such as mesoscopic capac- 
itors [3] and quantum dots [4]. The delay times took 
more attention in chaotic systems and random poten- 
tials to provide a statistical description of the time spent 
by wavepakets in the interacting region and the fluctu- 
ation of physical observables directly related to, such as 
thermopower [•">]. The distribution of the delay time was 
obtained in the case of one-dimensional disordered sys- 
tems by means of invariant embedding formalism in [6] 
and also in [7]. The more general result including the 
different classes of symmetry is obtained in [3] but still 
for ID systems (one mode). The next step of treating 
the N-modes scattering problems, is more complicated. 
However, the authors of [8] succeeded and gave the dis- 
tribution of individual partial delay times defined as the 
energy derivative of the phase shifts T P artlal = <^Ll. The 
joint probability distribution of the proper delay times 
was finally obtained in [9] which, expressed in terms of 
the inverse proper delay times 7i = — , reads: 



(1) 



/3 stands for the symmetry class (f3 = 1, under time re- 
versal symmetry (TRS), (3 = 2 broken TRS, and /3 = 4 
broken spin- rotation symmetry ). A is the mean level 
spacing. We need to mention that this result is obtained 
for large Hamiltonians(A r — > oo) and at Fermi energy ly- 
ing in the bulk spectrum(A ~ l/iV)[10]. In this paper, 
we give the distribution of Q at the edge of the spectrum, 
for small values of the Hamiltonian size. To do this, we 
start with the following representation of the scattering 
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FIG. 1. Tight binding model for a chaotic cavity with N = 7 
sites and M = 4 leads. In order to obtain Scattering matrix 
from circular ensembles CE, N should be larger or equal to 
M. The case N = M is referred here as a minimal chaotic 
cavity. 



matrix: 



[I - iirW^RW] x[I+ iirW^RWY 



(2) 



where we introduced the reaction matrix R = (E — A — 
R)~ l . R. is an N x N Hamiltonian matrix describing the 
closed chaotic system and W is M x N coupling matrix 
between the cavity and the conducting modes. This way, 
the matrix X = A — inWW^ represents the self-energy 
of the leads (conducting modes) where the real part A 
and the energy dependence of S will not be neglected 
as in the wide band limit WBL which is by the way the 
case of Eq. (1) obtained by [9]. Hereafter, the case of 
uniform semi-infinite leads is assumed, and this makes, 
with the independence and equivalence assumption be- 
tween modes, £ = (E/2 - - (E/2) 2 ) x I M [U]- We 
work also within the equal a priori probability ansatz by 
requiring uniform distribution for S which is the defini- 
tion of circular ensembles CE where (S) = 0. Cavities 
with this condition will be called completely chaotic. It 
is easy to show [13] that if the Hamiltonian is Lorentzian 
with center E — A and width T/2 — — we get iS cir- 
cular with the corresponding class of symmetry (fi = 1,2 
and 4) [5]. It is important to mention that this is true 
for all N > M so that the size of R. is not required to 
be infinite and thus, all the cavities with different sizes 
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are equivalent and lead to the same distribution of any 
observable expressed using the elements of S such as the 
conductance [12] or shot noise. Unfortunately, this is not 
true for observables containing energy derivative of S[ ] 
so that more attention is needed. The case N = M still 
ensures S from CE and we will call this case as minimal 
chaotic cavity. This case is very interesting since all the 
matrices become square and iS is diagonalized with the 
same energy independent rotation that diagonalizes T-L. 
Moreover, it allows to express || as a function of S. All 
in all, we obtain a simple expression for Q (N = M): 



Q 



r 



■Q~ 



i 



1 

f ' 



(3) 



where we put Q+ = and Q~ = §^f-. T' is the 

energy derivative of T. Before we give the distribution of 
Q, we need to start with the joint probability distribution 
j.p.d of Q + eigenvalues {r^}. We use a parametrization 
of type Eq. (2) and define: 
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and d = f3(N — l)/2 + 1. Here, the {e^} are the eigen- 
values of the matrix ttW^RW which is also distributed 
according to Lorentzian ensembles [5] [13]. After straight- 
forward calculations based on the properties of the Dirac 
delta function, we obtain: 
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J2 n»<j \Sif(Tj,Tj) - srf{T h Tj)f ^ 
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where f(x,y) = y(l — x)(l + y) and the sum is taken 
over all the ensembles of Sk taking values ±1. Proceed- 
ing the same way, we prove that Q~ exhibits the same 
distribution. We remind that since Q, Q + , Q~ , S and T-L 
are diagonalized by the same rotation, the distribution of 
the eigenvectors is the same as that of the eigenvectors of 
"H, ie uniform. That is why we restrict the study to the 
j.p.d of the eigenvalues. Going back to Q, we notice that 
at the half filling limit E = (the conduction band with 

uniform leads is [— 2to, 2to][15]), we have Q = and 
therefore, we immediately obtain from Eq.(5) the j.p.d of 
the eigenvalues {t^} of Q at E = just by changing in 
Eq. (5) / by / defined as f(x,y) = \/ (1 — x)y. We need 
now to generalize this result to any arbitrary energy. For 
this task, we rewrite the matrix Q at arbitrary energy as 
follows: 
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a = (f) 2 an d cos9 = l/\/l + T' 2 , T' being the energy 
derivative of T. The expression of the distribution of the 
matrix Q comes straightforwardly and can be put in the 



following form: 

™^/«^^-(^)>^ 

(7) 

where we defined here S = e + S. We used also in writ- 
ing Eq.(7) the invariancc properties of the Haar measure 
under the transformation S — > VSV' where the arbitrary 
unitary matrix V satisfies V = V T in presence of time 
reversal symmetry (3 = 1 and V = V R for broken spin- 
rotation symmetry j3 = 4. The letters T, and R denote 
respectively the transpose and the dual operations. One 
can notice in Eq.(7) the presence of the expression of Q at 

the half filling limit (Q = x ~3 ) this allows us to express 
the distribution Ve(,Q) at arbitrary energy as a function 
of the VE^oiQ)- O nc obtains from Eq.(7) the following 
formula: 



V E (Q) « \v E 



Q-1/T + a/2 



(8) 



This means that all what happens at arbitrary energy 
in the conduction band can be obtained from the for- 
mula obtained at the half filling limit. The j.p.d of the 
eigenvalues {f%} at arbitrary energy reads: 

n i<7 - \Sifa,b( T i,Tj) ~ Sjf a , b {Tj,Ti)\ P 



P{{n}) « 



E 

{...s k =±} 



Ui fa,b(Ti,Ti) 



(9) 

The function f a ^ is defined as: fa,b(x,y) = 
y/(b — x)(y — a). As usually for an open system, the pa- 
rameters a and b are directly related to the conducting 
leads by the broadening T: a = ^ 



and b 



This result worths to stop on it to make some comments. 
First, we find that the repulsion term \n — ry| is absent 
in a sense that when = Tj the j.p.d of the proper delay 
times does not vanish as observed in Eq.(l). This can 
be understood if one notices that for this case (N = M) 
the transformation between the eigenenergies {Ei} and 
the proper delay times {r^} is not one-to-one. So, even 
though it is true that having two energies Ei = Ej is for- 
bidden by the symmetries corresponding to = l,/3 = 2 
and j3 = 4, the situation Ei = —Ej is allowed, and this 
leads to r, = Tj . Second, we notice that the delay times 
are bounded and lie in the interval of length a which be- 
comes larger and close to infinity at the bottom of the 
conduction band (the continuum limit). We can notice 
also that the proper delay times, as defined here, can 
be negative and the only situation where they are al- 
ways positive is at the half filling limit where the WBL 
is equivalent. This causes no problem since we know 
that a potential can induce a delay or an advance of the 
wavepacket[2]. 

The numerical simulations of the case of two modes are 
shown in Fig. 2. The procedure to obtain these figures 
consists in sampling Hamiltonians from Lorentzian en- 
sembles, expressing the scattering matrices using Eq. (2) 
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FIG. 2. Joint probability distribution of the proper delay times for a minimal chaotic cavity(M = 2, N = 2) for different 
symmetries: left (9 = 1, middle /3 = 2 and right /3 = 4. All the scattering matrices belong to CE. The numerical data are 
obtained by sampling Hamiltonians from Lorentzian ensembles. The analytical result is plotted as well and the matching is 
perfect. Here the Fermi energy Ef = — lto. 



and then calculating the corresponding Wigner-Smith 
matrix Q. A histogram of its eigenvalues is done at the 
end. The analytical result P(ti,t 2 ) of Eq.(9) is plotted 
at the same time to check if the plotted surface passes 
from all the points obtained by sampling. The matching 
between numerical and analytical result is excellent. 
The Hamiltonian H describing the chaotic cavity with 
a Lorentzian distribution exhibits a density of states 
given by p(E) = f (r/2) . + g! gF+A) , [13] [16]. This for- 
mula simplifies too much at the Fermi energy, where 
the electronic quantum transport is pertinent, to become 
p(Ep) = which means that the mean level spacing 
at the Fermi energy is A = j^. Eq.(l) was obtained 
for a situation corresponding to T/2 = 1 and N — > oo, 
that is to say, for a very small mean level spacing cor- 
responding to the universality of the bulk spectrum[10]. 
In this universality class, the use of the Lorentzian en- 
sembles is not pertinent and it can be replaced by the 
equivalent Gaussian ensembles [13] offering other proper- 
ties such as the independence of the Hamiltonian ele- 
ments. If we consider now a mesoscopic system with a 
finite number of degrees of freedom (N finite) [14], and 
work with a Fermi energy lying in the bottom of the con- 
duction band (Ep ~ —2to), the density of eigenvalues 
p{Ep) tends to vanish whereas the mean level spacing A 
diverges (because r — > 0). This situation corresponds to 
the edge of the spectrum with a physics different from 
that of the bulk. We want to show that the system we 
considered in this paper (N = M) and we called it min- 
imal cavity is very important to describe this physics. 
Suppose now that the Hamiltonian % of the chaotic cav- 
ity is of size N x N and that the number of conduct- 
ing modes is M < N. We can not apply the method 
and the results obtained so far because % and iS are 
not diagonalizable by the same set of eigenvectors any- 
more and moreover the matrix diagonalizing S is en- 
ergy dependent and thus we do not obtain Eq.(3). To 
study such systems we prefer to use a decimation method 
which consists in taking out the whole sites describing 
the chaotic cavity and renormalize the sites (and links 
between them) to which they were attached (black sites 
in Fig. (1)) [17]. This is translated by taking a new effec- 
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FIG. 3. Density of eigentimes of chaotic cavities with N sites 
and M leads at the spectrum edge. The scattering matrices 
belong to circular ensembles CE with the different symme- 
tries ((8 = 1, /3 = 2 and /3 = 4). ~k corresponds to presence of 
time reversal symmetry (TR) ie j3 = 1 with M = 2, N = 23). 
♦ is obtained in absence of (TRS) ie /3 = 2 (with M = 9, 
N = 15). • concerns broken spin rotational symmetry ie 
£ = 4 with (M = 5, N = 10). The results are obtained by 
sampling ~H from Lorentzian ensembles with the correspond- 
ing /3-symmetry. Fermi energy Ef = — 1.999io- 



tive Hamiltonian %{E) = t e ^^, t^ with a M x N matrix 
Tij = Sij,i < M. This procedure helps in lowering the 
dimension of the matrices we work with to M x M size 
and thus the self energy S becomes proportional to iden- 
tity. This procedure is exact (up to non pertinent phase 
factor) and ensures the same Lorentzian distribution for 
H with the same center £ = Ep — A and width T/2 
[13] [5] which means that the final S still belongs to CE. 
Replacing % by H(E) in Eq.(2) still do not allow to re- 
cover Eq.(3) because the new effective Hamiltonian %{E) 
is energy dependent. What can someone notice is that 
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the derivative 



E-H 



t* becomes much 



P = 2 



more simpler at the bottom of the conduction band: In 
fact, wc get r — > which means a very narrow distri- 
bution, almost deterministic, and this allows to replace 
H by its mean £ in the derivative of H(E). Moreover, 
l/(E—£) = 1/A — > —ImxM in the vicinity of the bottom 
of the conduction band. This result entails the following 
important conclusion: whatever is the finite size N x N 
of TL, we end with S and ^| independent of N and there- 
fore, the result obtained with the minimal chaotic cavities 
(N = M) applies for all of the other situations(iV > M). 
Moreover , this goes beyond the time delay matrix, to 
all the observablcs which can be expressed using the el- 
ements of S or ^ [5]. We stress that, the larger is N 
the more T needs to get closer to zero. To illustrate this 
theory, we can look at the mean density of proper delay 
times for the model with N = M, and compare it to 
the case N > M at the edge of the spectrum. This mean 
density is defined as p(r) = Y^i(^( T ~ T i))- After straight- 
forward calculations we obtain the following result (see 
Appendix A): 

M 

fa) = m f ( \ (10) 

The result is /3-independent which goes back to the fact 
that the density of eigenvalues of the Lorentzian Hamil- 
tonian is /3-independent too[13]. The numerical simula- 
tions plotted in Fig. (3) confirm this result and show that 
at the edge of the spectrum, the formula also describes 
the case TV > M as it was expected. 

The Wigner time and the two modes scattering 
problem: The Wigner time t w is defined as the trace 
of the Wigner-Smith time delay matrix: t w = Tr(Q). 
This observable is very important since it can be di- 
rectly related to the change in the local density of states 
[18] [19] and thus enters in the evaluation of the ther- 
modynamic expectation of physical quantities where we 
take into account the effect of the coulomb interactions 
in the grand canonical average [4]. For the minimal 
chaotic cavity (N = M), the use of the n- level cor- 
relation functions of circular ensembles [20] [21] gives ac- 
cess to the first cumulant which is found (3- independent, 
(t w ) = ^(see appendix B) as well as the second cumu- 
lant (St*) = ^(|i| 2 )(see Appendix C) where t = Tr(5). 
For /3 = 1 we get(|i| 2 ) = -j^j and for (3 = 2 we have 
(\t\ 2 ) — 1 [21]. The whole distribution is more complex 
to obtain and can be expressed only in the simple cases 
N = 1 and N = 2. For N = 1 the result is immedi- 
ate: Virtu) = — , 1 The case N = 2 is more 

tt-^/ (b-T w ){T w -a) 

complicated but someone can notice that in this case the 
Wigner time is closely related to the Seebeck coefficient 
Sk [22] so that by using the integrals and techniques ex- 
ploited in [5] [22] we find the distribution at arbitrary en- 
ergy (M = N = 2): 
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FIG. 4. Distribution of the Wigner time in the absence of 
time-reversal symmetry (,3 = 2). Green plot is obtained nu- 
merically whereas the red curve is obtained with Eq. (11). 
The Fermi energy (here Ef — 0) lies in the edge of the 
Hamiltonian spectrum (p(Ef) = 2.1CF 4 ). We used the new 
self-energy E' to reach this density at Ef = 0. The size of 
the cavity is N = 26 with two conducting modes (M = 2). 



2-F1 is the Gauss hypcrgeometric function[2-'!]. This 
distribution is singular and has a logarithmic diver- 
gence near the mean Wigner time : Ve(t w ) ~ 
Cp log( ^ Tw ~^ Tw ^ ), where Cp is a /3-dependent constant. 
As for all the observables depending on S and |g , at the 
spectrum edge when T — > 0, Eq.(ll) describes very well 
the Wigner time distribution even for N > M = 2 (but 
N finite). 

The form of the self energy of a uniform lead implies 
that the edge of the spectrum is obtained for a Fermi 
energy lying in the vicinity of the bottom of the conduc- 
tion band. In order, to check that the physics discussed 
here is due to the spectrum edge (and not to the bot- 
tom of the band) we suppose that in addition to the uni- 
form semi-infinite lead there are stumps[24] attached to 
the cavity. This implies a new self energy £' = q + S. 
The additional part q, due to the stump, can be cho- 
sen energy-independent in some limits of the coupling 
term. When q ^> T/2, the level density at Ep becomes 
p(Ep) ~ ^ 1- This very low density implies that 
at any Fermi energy, we are in the tail of the spectrum. 
Since the energy derivative of £' does not change from 
the previous case, all the conclusions claimed so far re- 
main unchanged and moreover they become true at any 
Fermi energy: The minimal chaotic cavity gives the dis- 
tribution of the transport coefficients (depending on S or 
|g of any cavity with a finite N at the spectrum edge. 
This is successfully tested with numerical simulations: 
Fig. 4 shows the comparison between simulation and the 
analytical result of the distribution of the Wigner time. 
The simulation is based on sampling Hamiltonians giving 
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rise to scattering matrices from the Circular unitary en- 
semble (/? = 2) according to the new self energy £'. The 
parameter q is chosen quite large to reach the edge of the 
Hamiltonian spectrum at Ep = 0. The distribution is 
symmetric around the the Wigner mean time f w = p. 

As a conclusion, we proposed a model of minimal 
chaotic cavities with scattering matrices from circular 
ensembles and gave the exact distribution of the proper 
delay times and the Wigner time. We showed that this 
model is pertinent to describe in general for N > M, the 
statistics of all the physical quantities depending on S 
or/and §^ at the edge of the Hamiltonian spectrum. 
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Appendix A: Density of proper delay times 

The density of proper delay times is defined as follows: 
p(r) = J2(S(r-r i )) (Al) 

i 

this means: 

. M 

i J k=l 

Using the parametrization leading to Eq.(4) and the 
equivalence between the proper delay times we can write: 

. m 
p(r) = M / S(t - f(ei))P({e k }) J[ de k (A3) 
J k=i 

f(x) = j^, {n = T j f ). f-\x) = 

The properties of the Dirac delta function allow to write: 

M 

= (1 _ r) 3/2 Tl/2 M J ^ ~ r\r))V{{e k }) H de k 



k=l 



PCE(f-\r)) 



(1 - tY^t 1 / 2 
(A4) 



where pce is the density of eigenvalues of Lorentzian en- 
sembles which is given by (in the half filling limit) [13]: 

, N M 1 
PCE(e) - —^—2 ( A5 ) 

From Eqs.(A4) and (A5), we obtain (at E F = 0): 

M 1 



P(r) 



7T v/r(l-r) 



(A6) 



We notice that the density of eigentimes do not depend 
on the symmetry class which goes back to the fact that 
the Lorentzian density of eigenvalues is /3-indepcndcnt. 
The generalization of this result to any arbitrary energy 
can be obtained following the same steps, or just by using 
Eq.(8). We finally get: 

P(r) = —^==== (A7) 
7T y/ (b-r r-a 



Appendix B: Wigner mean time 

We show here how we calculate the Wigner mean time. 
We will use the fact that the marginal distribution of the 
eigenphases (phases of S eigenvalues) is uniform. We 
have: 

„ M 

(r w ) = / (n + t 2 ■ ■ • + T M )V({n}) J] dn (bi) 

. M 

M nV({n}) Y[dn (B2) 
J i=i 

. M 

M I Tifa)H{<Pi})n*Pi ( B3 ) 

l—l 

mJ^tM^ (B4) 



because J 'P({(p l })Yl i _ il dip l = ^[20]. We recall that 
for circular ensembles, we have V({ifi}) oc Y 



%< 3 i 

From Eq.(6) we deduce the relation between Tj and ipi'. 

a e +i{v>i+8) + e -K<Pi+e) 1 



n ] 



2 r 

a 1 
cos(<^ +0) + - (B5) 



2 ' r 

Eqs. (B5) and (B4) lead to the final result: 

, , M 



' W — \ 'W I 



r 



(B6) 



Appendix C: The second cumulant 

We would like to calculate the second cumulant of the 
Wigner time defined as : 

{5tI) = {{t w -t w Y) (CI) 
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It is better to do the calculation using the eigenphases 
representation as in Eq.(B5) : 



(C2) 



First, we notice that a translation ipi + 6 — > ipi does not 
change the result. Second, the expansion of this expres- 
sion with the use of the symmetries between eigenphases 
allow to write (i, j are arbitrary with i ^ j): 



{5t 2 w ) = T (iV((cosV)) +N(N- l)((cos(^)cos(^)) 

(C3) 

The first term is simple since the marginal distribution 
of one eigenphase is uniform [20]: 

M 



Rt(v) = j8(<p-<pi)V({<p k })f[d<p k = ± (C4) 
this gives : 



k=l 



((cos 2 (^)) 
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The second term involves the density-density correlation 
function. We have: 

. M 
R 2 (<p, tp') = S(tp- tpi)6{ip' - ipj)V({ip k }) Y[ dtp k 
J fc=i 



(C6) 



with (i 7^ j). 



Rsfay') = (6(<p- tpi)8(tp' - <pj)) 

= n(n-i) Si #j(^(y - f i)&W - <Pj)) (C7) 

and p((p) = jj Y^i 8(<f — fi)- We defined A = (A). 

One can notice that p(cp) = Ri(ip). We can express p{tp) 

otherwise 

P(<P) = w'E,?=l S (<P-<Pi) 

= ij2 N J2 n=+oc e m ^-^) (C8) 

2-kN Z-^n— — oo ^ ' 

where we put t n = Tr{S^ n ) = e" m ^ . The density- 
density correlation function is : 



p(<p)p(<P') 
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\27riV 



E^ 



|2pi»(ipV) 
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We want to express the following integral: 



s(tp) cos = J cos(ip) cos((p')R2((p, (p')d(fd(p' 

(CIO) 



This will be done term by term: 
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2^v) E^™' 2 / cos(^) cos(^')' 



in((p — <p') 



( l 



V 2ttN 



2ttN 
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The second term to evaluate leads to 



cos((/?) cos(y>')<5(y> — ip')p{ip)dtpdip 



cos (ip)p(cp)dtp 



(C12) 

All in all, we obtain the following result correct for the 
three ensembles(/3 = 1, (3 = 2 and /? = 4) : 



{8rl) = ^\t\\, 



(C13) 



for p = 1 |ti| 2 = 1 and for (3 = 2 the form factor is |ti| 2 
j^TTf- In the case /3 = 4, the form factor |£i| 2 = 2 n-1 ' 
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